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ABSTRACT 


A small eccentricity analysis was performed for a bearing with two feeding planes, 
each of which is assumed to be a line source. Numerical results were obtained for a 
range of bearing number, pressure ratio, feeding parameter, and orifice recess volume. 
A digital computer program was written to obtain these results and is included. Steady- 
state load and attitude angle were obtained, as well as stability data. Stability decreased 
markedly with increasing recess volume; moreover, for large recess volume and low 
bearing number, an increase in pressure ratio decreased stability. There was no corre- 
lation between stability and steady-state attitude angle for any of the cases studied. Fair 
agreement was obtained with available experimental data. 



STABILITY ANALYSIS FOR UNLOADED EXTERNALLY PRESSURIZED 
GAS-LUBRICATED BEARINGS WITH JOURNAL ROTATION 
by David P. Fleming, Robert E. Cunningham, and William J. Anderson 

Lewis Research Center 

SUMMARY 

A small eccentricity analysis was performed to predict the stability of an externally 
pressurized gas-lubricated journal bearing with two feeding planes, each of which is as- 
sumed to be a line source. Numerical results were obtained for a representative range 
of bearing number, pressure ratio, feeding parameter, and orifice recess volume. These 
results were obtained by using a digital computer program which is presented in an ap- 
pendix. Steady-state load and attitude angle were obtained, as well as stability data. 
Stability decreased markedly with increasing recess volume; moreover, for large recess 
volume and low bearing number, an increase in pressure ratio decreased stability. There 
was no correlation between stability and steady-state attitude angle for any of the cases 
studied. Fair agreement was obtained with available experimental data. 


INTRODUCTION 

Self-excited instability is one of the severest problems encountered in operating gas- 
lubricated journal bearings. At high speeds and low loads, the journal may precess about 
its steady-state position at a frequency which, for a self-acting bearing, is usually about 
one-half of the frequency of rotation. Hence, this phenomenon has become known as half- 
frequency whirl. In bearings other than the plain self-acting bearing, particularly in the 
externally pressurized bearing, the whirl frequency can be considerably different than 
one-half of the rotational frequency. Therefore, the phenomenon is more properly called 
fractional frequency whirl. Since fractional frequency whirl is excited by forces gener- 
ated within the fluid film of the bearing, it is not possible to eliminate it by increasing 
speed, as is the case with a synchronous resonant condition. Any attempt to increase 
speed once fractional frequency whirl has begun usually results in rubbing contact be- 
tween the journal and the bearing. 



A major part of the research in gas -lubricated bearings has been directed toward de- 
velopment of bearing configurations that will operate stably when the load is small or 
zero. Three designs have received primary consideration: herringbone -grooved, tilting- 
pad, and externally pressurized bearings. 

(1) Herringbone -grooved bearing. Shallow helical grooves are cut in either the jour- 
nal or the bearing to create an inward-pumping viscous pump. This type of bearing shows 
good promise for operating stably to high speeds (ref. 1), but clearances must be kept 
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very low, approximately 0. 5x10 centimeter per centimeter of shaft radius. These low 
clearances may cause problems at assembly or when thermal distortion is present, and 
they require extreme cleanliness of the lubricating gas. 

(2) Tilting -pad bearing. This is perhaps the most stable type of gas bearing (ref. 2), 
but the flexures or pivots required for each of the pads increase the complexity of the sys- 
tem and can lead to fatigue or fretting problems in long-term use. 

(3) Externally pressurized bearing. This type of bearing operates well at moderate 
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clearances (10 cm/cm) and has fairly good stability characteristics. The chief disad- 
vantage is the need for a source of lubricant under pressure. However, sufficient pres- 
sure may be generated within the system which uses the gas bearings, so that some of the 
high-pressure gas may be bled off and used for the bearings without imposing severe pen- 
alties on cycle efficiency. This is especially true in large systems. In these systems, 
externally pressurized bearings may be advantageous and even mandatory, because self- 
acting bearings large enough to support the required load are impractical. This is so be- 
cause the load capacity of a self-acting bearing is proportional to the square of its diam- 
eter, while the load imposed is a function of the rotating mass, which varies as the cube 
of the rotor diameter. 

Startup problems are common to all of these bearing types. Proper selection of 
bearing materials will minimize these problems. If the rotor can be "floated” by exter- 
nal pressurization before rotation is attempted, however, the materials problems is re- 
duced considerably. This is easily accomplished in the externally pressurized bearing, 
and it is common practice, even with tilting-pad bearings to include provision for exter- 
nal pressurization at startup (ref. 3). Obviously, the externally pressurized bearing is 
simpler than the tilting-pad bearing when the latter must be externally pressurized for 
startup. 

Since externally pressurized bearings may be the best choice for some applications, 
it is important to know the conditions under which these bearings will operate stably. Be- 
cause of the large number of parameters that may vary in a physical system, it is slow 
and expensive to rely solely on experimentation to gain needed information, and it be- 
comes extremely desirable to apply analytical techniques to the stability problem. 

A previous analysis (ref. 4) predicts stable regions of operation for externally pres- 
surized gas -lubricated bearings operating at finite (nonzero) eccentricities. Account is 
taken of orifice recess volume, and an attempt is made to account for having a finite num- 
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ber of orifices rather than a line source, which is usually assumed. However, the method 
is quite cumbersome and expensive to apply in terms of computer time. In addition, the 
analysis requires modification to be applicable to bearings in which the lubricant is intro- 
duced at more than one axial location. 

A need exists, therefore, for a technique that will provide stability information for an 
unloaded bearing. Pan (ref. 5) has outlined a general method that is applicable to both 
loaded and unloaded bearings . For the case of an unloaded bearing, the method is quite 
simple to apply: the only information needed is steady-state operating data for small am- 
plitude constant whirling. On the other hand, Lund (ref. 6) has derived the expressions 
for the steady-state load and attitude angle for an externally pressurized bearing operating 
with a small eccentricity. 

The analysis of reference 6 may easily be modified to yield a solution for steady 
whirling, including the effect of orifice recess volume. When this modified analysis is 
subjected to the stability criteria of reference 5, the regions of stable operation of an ex- 
ternally pressurized bearing are quite easily determined. 

The objectives of this report are to (1) combine the methods of references 4 to 6 to 
obtain a stability analysis for an externally pressurized gas -lubricated journal bearing 
operating at zero steady-state eccentricity, (2) develop and present a computer program 
to carry out the analysis, (3) present representative results of the analysis, and (4) com- 
pare the analytical results with experimental data. 


ANALYSIS 

The analysis is accomplished by applying the expressions from reference 4 concern- 
ing orifice recess volume to the analysis of reference 6 to obtain solutions for steady cir- 
cular whirling as in reference 7. These solutions are then subjected to the stability cri- 
teria of reference 5. The analysis is presented in its entirety for better unification of the 
various constituents. 

The configuration to be analyzed is that of an externally pressurized gas-lubricated 
journal bearing that has two circumferential rows of orifices (fig. 1). Although, in prac- 
tice, there are a finite number of orifices in each of the two rows, in order to make the 
problem tractable it is assumed that each row can be represented by a line source. The 
starting point in the analysis is the time -dependent, compressible, isothermal Reynolds 
equation (ref. 8): 
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All symbols are defined in appendix A. This equation may be made dimensionless by 
using the variables 


e = C = 


R 


R 



( 2 ) 
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Noting that P(3P/3x) = (1/2)(3P /3x), etc., the dimensionless form of equation (1) is 


A(h 3 -L(h 3 — V 2A — (PH) + 2o — (PH) (3) 

30 \ 30/ 3? \ 3? / 30 3 r 


in which the bearing number A is defined by 


A = 6 W R f 

P a \C) 


(4) 


and the frequency number a by 


a 
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The following two cases will be considered: 

(1) There is no whirling (co p = cr =0), but an externally applied load causes the bear- 
ing to operate with a small eccentricity ratio e . 

(2) There is no external load, but the journal is whirling in a small circular orbit 
about the bearing center with frequency co p . 

For both cases there will be a force generated in the fluid film which will tend to move the 
journal center toward the bearing center. 

If the journal and bearing axes are assumed to be always parallel, the dimensionless 
film thickness H for case 1 is given by 


H(0) = 1 + e cos 0 


and for case 2 by 


H(0, r) = 1 + e cos (9 - r) = 1 + e cos 9* 


As in reference 7, a rotating coordinate system is introduced by 

* 

9=9-r 


Then 


9 _ 9 
96) 90* 


and 
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For both cases 1 and 2, 


the Reynolds equation becomes 
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where, in the last term, 1 + e cos 0* has been substituted for H. 

The pressure P may now be represented as a power series in the eccentricity 
ratio e. For small values of e, it will be sufficient to consider only the first two terms 
of the series: 
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( 11 ) 


P(e*,S) = p 0 (O +eP 1 (e*,0 

If equations (11) and (7) are substituted into equation (10) and the resulting expression is 
considered an identity in e, a separate equation may be written for each power of e 
which appears. If powers of e higher than 1 are neglected, two equations result: 


32p 0 
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Pq sin 0 


= 0 


(13) 


Before complete solutions of these equations can be obtained, boundary conditions for 
Pq and Pj must be found. 

At the end of the bearing, the film pressure p equals the atmospheric pressure p . 
The conditions on Pq and Pj are then 

Pq = 1 and Pj = 0 at £ = 0 (14) 

The pressure is symmetric about the bearing midplane; thus, 

0P n 0P t 

— -0 = — - = 0 at C = (15) 

0£ 0? 


where 


L - L n 

k = 

D 


and 


«1 



(16a) 


(16b) 


With this condition, it is necessary to obtain a solution in only one-half of the bearing. At 
the feeding planes, the gas flow through the orifice restrictors equals the flow out through 
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the bearing. The dimensionless mass flow through a single orifice is given by 



To approximately account for the effect of inherent compensation, Lund (ref. 4) proposed 

r 2 -|i/ 2 

to modify the mass flow in equation (17) by dividing by the factor 11 + (5/H) J , where 
6 = a^/dC is the "inherent compensation factor. " 

The orifice recess (fig. 1), which has a volume V, can act as a reservoir for the 
supply gas . The actual mass flow into the bearing from one orifice recess is 


P a V c cu D 0P_ 
= M - a — - P — £ 
D c 3 t 


(18) 


in which P c is the dimensionless pressure in the recess. Since the inherent compensa- 
tion factor is assumed to be small, the pressure in the bearing adjacent to the orifice re- 
cess will be infinitesimally different from the pressure in the recess, and is also denoted 
by P . The dimensional orifice mass flow modified for inherent compensation is denoted 
by M c and is given by 


M c 


Ja2p sPa 

{a t ) 1 / 2 



(19) 


The Hagen -Poiseulle law gives the flow in the bearing away from the feeding plane. This 
must equal the flow out of the orifice recesses. Since there are N/2 orifices in each 
row, equations of mass flow into and out of the bearing yield, per unit of circumference, 
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( 20 ) 
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The right side of this expression denotes a discontinuity in the axial derivative of the pres- 
sure at £ = £ due to the gas entering the bearing at this point. The mass flow may be ex- 
pressed in a power series in e in the same way as the pressure. If, again, only powers 
of e through 1 are retained, 


9M„ 


M c = M cQ + 6 


9e 


e=0 


( 21 ) 


If one now substitutes the expression found by differentiating equation (19) for 9M c /9e in 
equation (21) and then combines equations (7), (11), and (21) with equation (20), an expres- 
sion results which, again, may be considered an identity in e. If powers of e only 
through 1 are retained, two equations result 


q = 


0± Q 

3? 


1?=!' 


3? 


k=H + 


( 22 ) 


3(P 0 P 1> 


. 8 < P 0 P 1> 

_ q [ 1 9m 

3 + 26 2 

3? 


3? 

2 \m n de 
\ 0 

e=0 1+62 


3( p 0 p i) 
cos e I + CM//J 

1 1 de* 




(23) 


in which the following definitions are used: 
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(24) 
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To evaluate 9m/9e, one writes 
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If \[/q is defined by 



(29) 


boundary condition (23) becomes 
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Now that the boundary conditions have been established, the differential equations 
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can be solved. The solution of equation (12) for Pq (the concentric solution) is easily ob- 
tained as 


Pq = i + q? 

Pq = 1 + 


C<! I 




(31) 


A solution of equation (13) may be obtained by assuming that 


P 0 p i = 


G(?)e 


- 10 * 


(32) 


in which G is a complex function of £. If this expression is inserted into equation (13) 
and it is noted that sin 9* = ), one obtains, after dividing through by e"^ , 


G" 


-G + i(A-o) ^ + P Q =0 


(33) 


The primes denote differentiation with respect to the argument The boundary condi- 
tions on Pj from equations (14), (15), and (30) become, respectively, 


G = 0 at £ = 0 

(34) 

= 0 at 5 = 1+ ^ 

(35) 


and 


- G* 


+ (i// Q + icn/^G 


4 


= q 


8* + 1.5 

i *2 
1 + 6 


(36) 


The ordinary differential equation (eq. (33)) may now be solved numerically by using 
a forward integration technique such as Runge-Kutta. The procedure is: Assume a func- 
tion G^ which satisfies all the conditions for G except that G'^(0) = j8, an arbitrary 
complex number. Integrate the differential equation (eq. 33)) for G^ from £ = 0 to 
£ = £ and apply equation (36). Continue the integration to £ = £ + ^ and store the value 
of G’^(4 + £j) = B. Now assume another function G y for which G’^(0) = y and repeat 
the above procedure to obtain G'^(4 + |j) = r. The linearity of the differential equation 
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(eq. 33)) allows the use of a linear combination of the two solutions just obtained to yield 
a third solution which satisfies boundary condition (35). Set G'(0) = (/3r - yB)/(r - B) 
and a third integration will yield the desired function G(£) . 


Determination of Load 


The radial and tangential forces acting on the bearing may now be determined by in- 
tegrating the film pressure over the bearing area: 



Since Pq is independent of 0 , its contribution to the bearing forces will be zero. Per- 
forming the rest of the 0 integration, noting that 


sfe * 

cos 0 + G^ sin 0 

(where G^ and G^^ are the real and imaginary parts of G, respectively), one has 
(in dimensionless variables) : 




The dimensionless forces in equation (37) are defined by 



(38) 


f = — 

r *P a LD 

► 


L <*a LD J 

Thus, the bearing load components are determined simply as the integral of G/Pq. This 
integration is conveniently carried out at the same time the differential equation (eq. (33)) 
is solved. In fact, fiG/P^dt, may be expressed as a linear combination of /(G^/P^d^ 
and f{ G^/Pq) 9? and only two integrations of differential equation (33) are necessary. In 
terms of G^ and G^, 



The bearing radial and tangential load components, which are illustrated in figure 2, are 



Figure 2. - Notation for eccentric bearing. 
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now determined as functions of bearing geometry, supply pressure ratio P , bearing 

O 

number A, and frequency number a. The total bearing load 


W = 



and attitude angle 


<P 


= tan 


-1 



( 40 ) 


(41) 


may now be calculated. 


Stability Criteria 

If the bearing is operating stably, o = 0. To determine the threshold of instability, 
according to the work of Pan (ref. 5), a is varied until F^. = 0. Under this condition, 
the centrifugal force due to the whirling mass is equated to the radial bearing force to find 
the neutral stability condition: 


Mew“ = F 
pn rn 


(42) 


The subscript n is used to denote the condition where F^. = 0. When the bearing is 
whirling, it is assumed that the flow of gas through the orifice does not vary with time, 
since the orifices lie behind the recess volume, which acts to dampen flow fluctuations. 
A dimensionless mass may be defined by 


M 


1 " 



(43) 


In terms of previously calculated quantities, for the neutral stability condition is 
given by 


M ln = 


36 f 


rn 


n 


(44) 
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It is shown in reference 5 that M ln is an upper limit of for stability if the quantity 
df^/da is negative at a - a n ; conversely, Mj n is a lower limit for stability if df^/do 
is positive at <7 = a . The definition of used herein differs by a factor of 2 from 
the dimensionless mass used in reference 4; however, it is identical to that used in sta- 
bility studies of herringbone-grooved gas bearings (refs. 1 and 9). 

Another convenient parameter expressing the mass in dimensionless form is 





MCco 2 

P a LD 


It is related to Mj by 


M 



(45) 


(46) 



(47) 


Each of these dimensionless masses, Mj and Mg, has its particular advantages. 
Examination of equation (43) shows that depends only on the properties of the bear- 
ing and lubricant; it is independent of speed. While Mg is a function of speed, it will be 
seen in the RESULTS AND DISCUSSION section that Mg for many bearing configurations 
is nearly constant for a wide range of A . 


RESULTS AND DISCUSSION 

Numerical results have been obtained for dimensionless load W, attitude angle cp, 
masses M^ n and Mg n , and whirl ratio Wp/w by using the procedure described in the 
preceding section. All computations were performed by a high-speed digital computer. 
The computer program, with sample input and output, is given in appendix B. The bear- 
ing configuration chosen has a length-to -diameter ratio L/D of 1.5, with two circumfer- 
ential feeding planes each located a distance L/4 from the center plane of the bearing; 

that is, L../D = 0.75. The inherent compensation factor 5 was set at 0.447, which makes 
, o 

1 + 6 = 1.2. These values correspond to those of an externally pressurized journal bear- 
ing being evaluated experimentally. In reference 10, orifice discharge coefficients are 
given as a function of the pressure ratio across the orifice and the Reynolds number of the 
flow. For the conditions encountered in externally pressurized bearings, a reasonable 
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mean value is 0.7. Accordingly, this value was used in the calculations. Three values 
of the feeding parameter A^. were used: 2, 4, and 10. 

Steady-State Results 

The steady-state results are identical to those which would be produced by the analy- 
sis of reference 6, except for the use of the inherent compensation factor? However, in 
that reference, no results are given for a bearing with double-plane admission. Inspec- 
tion of differential equation (33) and boundary conditions (34) to (36) together with the 
definition of ^ (eq. (26)), shows that the recess volume ratio v enters in only through 
\p^, which appears in equation (36). For steady-state operation, a , which multiplies \f/^, 
is zero, and thus the solution does not depend on v. Accordingly, plots of dimensionless 
load and attitude angle, in figures 3 and 4, are valid for all values of recess volume ratio, 



Figure 3. - Steady-state load capacity of externally pres- 
surized bearing. Ratio of bearing length to diameter, 

1. 5; ratio of distance between feeding planes to bearing 
diameter, 0.75; orifice discharge coefficient, 0.7; in- 
herent compensation factor, 0.447. 
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Figure 4. - Steady-state attitude angle for externally pres- 
surized bearing. Ratio of bearing length to diameter, 

1. 5; ratio of distance between feeding planes to bearing 
diameter, 0.75; orifice discharge coefficient, 0.7; in- 
herent compensation factor, 0.447. 


Figure 3 shows the load capacity increase due to rotational effects. This increase is 

larger for larger values of the feeding parameter A^., but is relatively insensitive to 

pressure ratio P . For higher values of P , however, the zero speed load is higher, so 
s s 

that the total load capacity at all speeds does increase significantly with P . For all but 

s 

the lowest value of P , the load capacity of the pressurized bearing exceeds that of the 

s 

unpressurized, plain bearing. For A^. = 10 and P g = 2, however, the unpressurized 
bearing carries a higher load than the pressurized bearing for 1/2 < A < 5. 

Attitude angles are shown in figure 4. The attitude angle cp at zero speed is always 
zero. Attitude angle increases rapidly with A initially, and then falls asymptotically 
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to zero as A — There are three curves in figure 4 which lie partly above that for the 
plain bearing, indicating that a pressurized bearing does not always have a smaller atti- 
tude angle than a self-acting bearing. An increase in pressure ratio always reduces cp, 
but the behavior with changing A^ is not so clearcut. For low pressure ratios and high 
A, (p increases with increasing A^., but the opposite trend appears to occur for low A 
and high P . 


Stability Results 

The dimensionless mass Mg n is plotted against bearing number A in figures 5(a), 
(b), and (c) for values of Aj. of 2, 4, and 10, respectively. In each of these figures are 
curves for P_ of 2, 5, and 10, and recess volume ratios v of 0, 0.2, and 0.4. The 
curves denote the largest value of M£ for which the bearing will be stable. 

The most striking feature of the figures is that Mg n does not vary with A when 
v = 0; that is, when there is no orifice recess, the stability as measured by M 2 n is in- 
dependent of the bearing number. For finite recess volumes, the dependence of 
on A becomes stronger as v increases; Mg n increases with A to an intermediate 
value of A, and then decreases. An increase in v always decreases stability. 

Another interesting result is that increasing the pressure ratio does not always in- 
crease stability, as is commonly assumed. At low values of A an increase in pressure 
ratio may actually decrease stability. This effect is accentuated more and more as the 
recess volume is increased. The reason is that high pressure ratio and large recess vol- 
ume are the conditions which promote pneumatic hammer. 

A comparison of figures 5(a), (b), and (c) shows that a variation in A^ has the same 
qualitative effect as a variation in P ; that is, an increase in A. increases stability at 
low values of v, and decreases stability at high values of v and low A. 

Some preliminary experimental data from the NASA Lewis externally pressurized 
test facility are plotted in figure 5(b) . The recess volume ratio in the experimental setup 
was approximately 0.4; other parameters are approximately as shown in figure 5(b). The 
pressure ratio is indicated next to each plotted point. The experimental points fall ap- 
proximately where curves would be for v = 0. 1, except for the point P = 2. 2, which is 

s 

lower, and, in fact, close to the line v = 0.4. Agreement can generally be considered 
fair. Part of the discrepancy between analysis and experiment could well be due to the 
number of orifices in the test apparatus. Whereas the analysis assumes a line source, 
there are, in fact, only six orifices in each row in the experimental setup. This could 
well reduce the effect of the large recess volume. 

The experimental apparatus contains two externally pressurized bearings with a rela- 
tively short span. Consequently, when the rotor became unstable, it whirled in a conical 
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Dimensionless rotor mass, M 


Pressure 

ratio, 

p s 



Orifice 

recess 

volume, 

v 


0 

□ ~ 4 (experimental) 


(a) Feeding parameter, A t , 2. 



<b) Feeding parameter, A t , 4. (c) Feeding parameter, A t , 10. 

Figure 5. - Stability limits of dimensionless rotor mass ZVU n for externally pressurized bearing. Ratio of bearing length to diameter, 1. 5; ratio of 
distance between feeding planes to bearing diameter, 0.75; orifice discharge coefficient, 0.7; inherent compensation factor, 0.447. 
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Dimensionless rotor mass, 


Pressure 



tance between feeding planes to bearing diameter, 0.75; orifice discharge co- 
efficient, 0.7; inherent compensation factor, 0.447; feeding parameter, 4. 



mode. The mass which was used in the experimental values of M 2n was an "equivalent 
mass for conical whirl" defined by 




) 


wherein I tr and I are the transverse and polar mass moments of inertia, respec- 
tively, and l is the bearing span. This form is suggested by the solutions to the equa- 
tions of motion of a rigid rotor presented in reference 11. 

The data of figure 5(b) have been replotted in figure 6 with M^ n as the ordinate in- 
stead of M 2n . As already mentioned, has the advantage that its definition does not 
include speed, so that in a given physical setup it remains constant, except for the varia- 
tion in clearance due to centrifugal force. The experimental points of figure 5(b) are also 
replotted, and the change in M-, due to centrifugal growth of the rotor is evident. 

The abscissa has been extended from A = 0. 5 to A = 0 on a linear scale. For a 
pressure ratio of 10 and recess volume ratio of 0.4, M^ n remains finite for A = 0, indi- 
cating a susceptibility to pneumatic hammer. For the lower pressure ratios and recess 
volume ratios, M^ n — °° as A — 0, indicating that the bearing is inherently stable when 
there is no rotation. The tendency toward pneumatic hammer may also be deduced from 
figure 5. If M 2n approaches a finite value as A — 0, there will be no hammer. On the 
other hand, if M 2n — 0 as A — 0, M^ n remains finite as A — 0, and pneumatic hammer 
is possible. 

The curves of figures 5 and 6 indicate the heretofore largely unrecognized importance 
of recess volume in journal bearings. Within the range of variables plotted, changes in v 
cause larger changes in stability limits than do pressure ratio or feeding parameter. In 
the design of bearings where instability may be a problem, therefore, it is essential that 
the recess volume ratio be minimized. 

According to the analysis, the instability will occur at some characteristic frequency. 
For finite A, this frequency may be conveniently expressed as a whirl ratio a> p /w. Fig- 
ure 7 illustrates the characteristic whirl ratios corresponding to the stability curves of 
figures 5(b) and 6. When the recess volume is zero, the whirl ratio is always 0. 5, or 
"half frequency. " Whirl ratio a> p /co increases with increasing v, and decreases with 
decreasing A, approaching co /w = 0.5 as A — °°. For the case of pneumatic hammer, 

Co = A = 0 and, of course, co p /w = °°. Thus, the values for P g = 10 and v = 0.4 in fig- 
ure 7, becoming large for small A, indicate again that hammer may be a problem for 
this configuration. 

Experimental values are also plotted in figure 7. Again, the pressure ratio is indi- 
cated next to each plotted point. These points bear the same relation to the analytical 
curves as in figures 5(b) and 6; that is, the experimental points correspond to analytical 
curves for smaller values of v than were measured. 
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Figure 7. - Frequency of instability. Ratio of bearing length to diameter, 1. 5; ratio of distance between feeding 
planes to bearing diameter, 0.75; orifice discharge coefficient, 0.7; inherent compensation factor, 0.447; feed- 
ing parameter, 4. 


Relation Between Attitude Angle and Stability 

It is popularly assumed that the steady-state attitude angle is a measure of the sta- 
bility of a fluid-film bearing, a low attitude angle indicating a stable bearing, and a high 
attitude angle a potentially unstable bearing (ref. 12). This is definitely not the case for 
the externally pressurized gas -lubricated bearing studied in this report, as an examina- 
tion of figures 4 to 6 will show. 

To eliminate the effect of variable recess volume, at first only the stability curves 
for v = 0 will be considered. Looking at the variation with A, one sees that, for all 
values of pressure ratio and feeding parameter, cp first increases with A and then de- 
creases asymptotically to zero, while M^ n decreases continuously and M 2 n remains 
constant. Thus, correlation between cp and stability is essentially nonexistent. Now, 
consider the variation with feeding parameter A^.. Stability increases with A^., while at 
low pressure ratios cp also increases with A t> At higher pressure ratios, the behavior 
of cp with A is anomalous. A plain (unpressurized) bearing, corresponding to A^ = 0, 
has an attitude angle lower than that for some pressurized bearings, except at very low 
A, and yet the plain bearing is always inherently unstable. When pressure ratio is con- 
sidered, the attitude angle decreases with increased P g while stability does increase. 
However, if now recess volumes other than zero are considered, stability does not al- 
ways increase with increasing P . One concludes that attitude angle is a very unreliable 
predictor of stability. 
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SUMMARY OF RESULTS 


A small eccentricity stability analysis has been performed for an externally pres- 
surized gas-lubricated bearing that has two feeding planes, each of which is considered 
to be a line source. A computer program was written to implement the analysis. The 
following results were obtained for a bearing having a length-to-diameter ratio of 1. 5 
with two feeding planes located midway between the center and ends of the bearing (quar- 
ter plane admission) : 

1. Stability decreases markedly with increasing recess volume. 

2. There is no correlation between stability and steady-state attitude angle. 

3. With large recess volumes and low bearing numbers, an increase in pressure ra- 
tio may decrease stability. For sufficiently high recess volumes and pressure ratios, 
pneumatic hammer is possible. 

4. The increase in steady-state load capacity due to rotation increases with increas- 
ing feeding parameter but is relatively insensitive to pressure ratio; however, total load 
capacity does increase with increasing pressure ratio. 

5. When the orifice recess volume is zero, the dimensionless mass Mj n does not 
vary with bearing number A . 

6. Steady-state attitude angles are sometimes larger than those in plain (unpres- 
surized) bearings. 

7. The whirl ratio at the threshold of instability is always 0. 5 when the recess vol- 
ume is zero. Whirl ratio increases with increasing recess volume, but approaches 0. 5 
asymptotically at large bearing numbers . 

8. Preliminary experimental data are in fair agreement with the analytical results. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 5, 1968, 

129-03-13-05-22. 
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APPENDIX A 


SYMBOLS 


a 

orifice radius 

P 

dimensionless pressure, p/p 

a 

C 

clearance at zero eccentricity 

P 1 

perturbed pressure 

D 

bearing diameter 

P 

pressure 

d 

orifice recess diameter 

Pa 

atmospheric pressure 

e 

eccentricity 

q 

A f m P 
t o s 

F 

bearing load component 

R 

bearing radius 

f 

dimensionless load component, 
F/ep a LD 

a T 

gas constant times absolute 
temperature 

G 

dimensionless complex function 

t 

time 


of ? 

v c 

orifice recess volume 

H 

dimensionless clearance, h/C 

V 

orifice recess volume ratio, 

h 

local film thickness 


NV c /ttDLC 

Vtr 

rotor polar and transverse mass 

W 

total bearing load 


moments of inertia 

X 

coordinate in direction of motion 

i 

yTi 

z 

axial coordinate 

k 

specific heat ratio 

a 

orifice (discharge) coefficient 

L 

bearing length 

P ,r 

complex values of G'(0) 

L 1 

distance between feeding planes 

b, r 

complex values of G'(i- + £j) 

l 

bearing span 

6 

inherent compensation factor, 

c\ 

M 

rotor mass per bearing 


a vdC 

M 

lubricant flow rate 

e 

eccentricity ratio, e/C 

Mi 

dimensionless rotor mass, 
(MP a /2L M 2 )(C/R) 5 

? 

dimensionless axial coordinate, 
z/R 

m 2 

dimensionless rotor mass, 

e 

angular coordinate 

m 

MCo> 2 /p o LD 

cL 

dimensionless lubricant flow 

e* 

angular coordinate, Q - o> t 

r 

bearing number, 6pa>(R/C) 2 ^p a 

N 

rate, M(#T) 1 / 2 y4ra 2 P g p a 
number of orifices per bearing 

A 
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A t 

feeding parameter, 

CO 

rotational speed 


LNa 2 («T) 1/2 ]/p C 3 (l + 5 2 ) 1/2 

CO 

p 

frequency of whirling 

M 

lubricant dynamic viscosity 

Subscripts: 


(L - L.)/D 

b 

condition immediately down- 

*1 

a 

X 

V D 

frequency number, 
12 M w p (R/C) 2 /p a 

dimensionless time, co p t 

c 

stream of orifice recess 

condition immediately down- 
stream of orifice 

T 

n 

0 

condition at which = 0 

zero eccentricity 

<P 

^0 

attitude angle 

-(A t /2P 0c ) [3m/3(P c /P s )]| £=o 

r 

s 

radial 

condition upstream of orifice 

*1 

nv c /.d 2 cp 0c 

t 

tangential 
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APPENDIX B 


COMPUTER PROGRAM TO DETERMINE STABILITY OF EXTERNALLY 
PRESSURIZED GAS-LUBRICATED JOURNAL BEARING 

The program for determining steady-state and stability characteristics of the ex- 
ternally pressurized bearing is written in FORTRAN IV, version 13, for use on the IBM 
7094-n digital computer at Lewis Research Center. Minor modifications may be neces- 
sary to allow the program to be used with other computing systems . 


Program Input 

Two variables are set within the program: k, the specific heat ratio of the lubricant 
gas, and RKEP, an accuracy parameter. As the program operates, the step size in the 
Runge-Kutta differential equation solver is successively halved until the magnitude of the 
bearing load varies by less than RKEP times W. 

The remainder of the information needed by the program is read from punched cards. 
One card contains the array of bearing numbers A to be used (up to 19 values) . A pair 
of cards shows each configuration to be analyzed. The first of this pair gives geometric 
and operating parameters; the second gives initial estimates of the frequency number a, 
one for each A value. Specific formats follow. 

Bearing number array card: Format (14, 19F4.0) 

First four-column field: the number of A’s in the array, integer format, right 
adjusted. 

Succeeding four-column fields: the values of A for which calculations are desired, 
real format (with decimal point) . 

Geometric configuration card: Format (10F8.0) 

Succeeding eight-column fields contain, in real format (with decimal point) : 

(L - Lj)/D, Lj/D, A t , P g , v, a, and 1 + 5 2 . 

Frequency number estimate card: Format (4X, 19F4.0) 

First four-column field is not read. Succeeding four-column fields contain, in real 
format (with decimal point) the initial estimate of frequency number a , one value 
for each A. 

Any number of pairs of geometric configuration and frequency number estimate cards may 
be used. 

Sample program input follows the program listing included at the end of this appen- 
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dix. The first card, the A array, indicates there are 19 A values ranging from 0 to 100. 
This card is followed by the first geometric configuration card, which indicates that 
(L - Lj)/D = 0.75, L x /D = 0.75, A t = 4, P g = 5, v = 0.4, or = 0.7, and 1 + 6 2 = 1.2. 
After this card is the first a estimate card showing initial estimates of a , which cor- 
respond to the A array, of 3 to 100. A second configuration card and a estimate card 
conclude the sample input. 


Program Output 

Output consists of two sheets for each configuration: A ’’working sheet” and a ’’cal- 
culation summary sheet.” 

The working sheet shows the geometric and operating parameters, trial values for 
determination of the pressure downstream of the orifices, final value of orifice down- 
stream pressure, i//q, and RKEP. For each A in the array, there appear, for each 
trial value of o : A, f r , f^, W/(ep a LD), W/je(P g - P a )LI>] , <p, Mg, final Runge-Kutta 
step size, M^, o>p/u>, and o . A sample working sheet appears as figure 8. The data ap- 
pearing thereon were generated from the values shown in the sample input. 

The calculation summary sheet again shows the geometric and operating parameters, 
and, for each A of the A array: W/(ep & LD), cp, Wp n /o>, c n , Mg n , and M^ n . Figure 9 
shows an example of a summary sheet. Again, the data correspond to those of the sam- 
ple input. Program execution time for the sample input was 0.9 minute. 
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$ I BFTC L PS 

DIMENSION SLAM I 20) ,SW(20> ,SPHI (20) , SW3 ( 20 ) , SCM ( 2 0 ) ,SSIG(20) 

DIMENSION ACM ( 2 0 ) 

COMMON .XI 

COMMON /CRK/ QD, XII* CG 
COMMON /CDR/ 0, GO 

COMMON /CFUN/ V, LT, ALFA, PCR, K 

COMMON /CHS/ LAMBDA, PSO , PS1 , DX , RDX ,CM A , V V , W , PHI , CMC, WR,RKEP 

REAL LAMBDA, LT, K 

COMPLEX DX, GD, CG 

DATA EPS/1. E-5/, EPSW/l.E-3/ 

K = 1.4 
RKEP = .005 

READ (5,1) IS, (SLAM(I), I = 1, IS) 

1 FORMAT (14,1 9F4 . 0 ) 

3 READ (5,5) X I , X 1 1 , LT , V , VC , ALF A , D1 

5 FORMAT (10F8.0) 

READ (5,2) ( SS I G ( I ) , I = 1, 1ST 

2 FORMAT { 4X , 1 9F4 . 0 ) 

VV = V 

WR I T E ( 6,4) 

4 FORMAT ( 1H1/8X, 6 7HP-PERTURBAT I ON ANALYSTS FOR HYBRID BEARING IN S 
1TEADY CIRCULAR WHIRL ) 

C X I =L ( OUTBOARD > / D , X I 1 =L ( INBOARD) /D, LT=L AMBDA (T> , ALF A=OR I F I CE COEFFICIENT 
C V=PS/PA, VC=FEEDER HOLE VOLUME RATIO, D1 = 1+DEL T A**2 
C XI1=0 DENOTES SINGLE PLANE FEEDING 

WRITE (6,6) XI,XI1,LT,D1,K,ALFA,V,VC 

6 FORMAT (6HL L/D 8X,4HL1/D 7X , 8HL AMBDA T 4X , 10H1 + DELT A**2 3X, 

1 7HK = CP/C V 7X , 4HAL FA 6X,7HV=PS/PA 7X,2HVC /11GI2.3) 

PCR = (2./(K+;.) >**<K/(K-1. ) ) 

YA = PCR 
WRITE (6,8) 

8 FORMAT (20HLS0LVE FOR YA = PC/V/) 

DO 20 I = 1,10 
CALL FUN(YA,DYA,PSO) 

WR I TE ( 6 , 10 ) I, YA, DYA 
10 FORMAT (14, 4X , 2 H Y A G14.6,3HDYA G14.6) 

IF( ABS (DYA/YA) .LE.EPS) GO TO 25 
YA = YA + DYA 

20 YA = AMIN1 (YA, .99) 

GO TO 3 

25 0 = ( (YA*V>**2 -1. )/XI 
OD = 0* ( D1 + • 5 ) / D1 
PSO = PSO/ ( 2 . ^YA’XV ) 

PS1 = VC*(XI+XI1 )/( YA*V> 

WRITE (6,26) YA, PSO, PS1 , RKEP 

26 FORMAT ( 9HK PC/PS G10.3,5H PSO G10.3,6H PS1 G10.3,7H RKDX 
1 25H REDUCED UNTIL ERROR .LE. G10.3) 

DO 300 1=1, IS 
WRITE (6,140) 

140 FORMAT 4 8HK LAMBDA 4X , 8HFR / EPALD 3X , 8HFT / E PALD 4X,7HW/EPALD 2X, 

1 1 2HW/E ( PS— PA ) LD 4X,3HPHI 4X,9HMCW2/PALD 7X,4HRKDX 5X, 

1 12HMPC5/2LR5MU2 IX, 11HWHIRL RATIO 1X,12HVIB NR SIGMA ) 

LAMBDA = SLAM(I) 

SIG1 = SS I G ( I ) 

DS - 2 . 
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IF (LAMBDA. GE. 15. ) DS=4. 

CALL HYST (0.,GI0) 

IF (LAMBDA. EO. 0. ) GIO=.25 
SW ( I ) = W 
S PH I ( I ) = PH I 
CALL HYST (SIG1, Gil) 

SIG = SIG1 

IF ( ABS(GI 1/GIO) .LE.EPSW) GO TO 299 
DS1 = SIGN (DS, G 1 1 ) 

SIG2 = SIG1 
DO 150 J = 1, 20 
SIG2 = SIG2 + DS1 
CALL HYST (SIG2, GI2) 

SIG = SIG2 

IF (ABS(GI2/GI0). LE.EPSW) GO TO 299 
DS2 - SIGN (DS, G 1 2 ) 

IF ( DS2 - DS1) 220, 210, 220 
210 SIG1 = SIG2 
150 GIL - GI2 
GO TO 299 

220 DO 235 J = 1 , 10 

SIG = (GI2*SIG1 - GI1*SIG2)/(GI2 - Gil) 

CALL HYST (SIG, GI) 

IF(ABS(GI/Gin). LE.EPSW) GO TO 299 
IF (ABS(GI2) .GT.ABS(GI1 ) ) GO TO 230 
GI1 = G 1 2 
SIG1 = SIG2 
230 G I 2 - GI 

235 SIG2 = SIG 

299 S W3 ( I ) =WR 
AC M ( I ) = CMA 
SCM( I )=CMC 
SS I G ( I ) = SIG 

300 CONTINUE 
WR I T E ( 6,4) 

WRITE (6,600) XI,XI1,LT, 

600 FORMAT ( 1HL 8X,3HL/0 8X,4HL1 
1 7HK=C P/C V 6X,4HALFA/5X,6G 
213HCL.VOL.RAT.VC 3X,5HPC/PS 
WRITE (6,325) ( SL AM ( I ) , SW { I 

1 1=1, IS) 

325 FORMAT ( / 1 HL 7X,6HBRG NR 3X,8HBRG LOAD 2X,22HATT ANGLE WHIRL RAT 
Ilf) 2 X , 6H V I B NR 3X , 20HST ABILITY PARAMETERS / 8X , 6HL AMBDA 4X, 

2 7HW/EPALD 5X,3HPHI 6X , 8HW3/0MEGA 4X,5HSIGMA 4X , 9HMCW2 /PALD 2X, 

3 12HMPC5/2LR5MU2 / ( / 7X , 2G 1 0. 3, G 11 . 3, G 12 . 3 , 2G 10. 3 , G 12 . 3 ) ) 

GO TO 3 

FND 

SIRFTC FUNF 

SUBROUTINE FUN (Y, DY, PSO) 

COMMON XI 

COMMON /C FUN/ V, LT, ALFA, PCR, K 
REAL L T , M , K 
IF(Y.GT.PCR) GO TO 10 

M = ALFA*S0RT(2.*K/(K+1. ) )*PCR**( l./K) 

PSO = 0. 

GO TO 15 

10 T = SQRT( l.-Y**( (K-l. )/K) ) 

TC = SORT (2.*K/(K-1. ) ) *AL FA 
M = TC^T^Y** (l./K) 


D1,K,ALFA,V,VC,YA,PS0» PS1 »RKEP 

/D 6 X , 8HL AMBDA T 4X, 10H1 + DELT A 8 ** 2 3X, 

12 .3//7X ,7HV=PS/PA 2X, 

8X , 3HP SO 9X , 3HPS 1 9X , 4HRKE P/ 5X , 6G1 2 . 3 ) 
) , SPH I ( I > , SW3 ( I ) , SS I G ( I) , ACM ( I ) ,SCM( I ) , 



(K+l.)/2.)/K 


PSO = -LT*TC/T*(Y**( (l.-K)/K) - 
15 1 - Y*V*V/XI 

F = LT*V*M - Y*T + l./XI 
OF = V*PSQ + 2.*T 
DY = F/DF 
RETURN 
ENO 

$ I B FTC HYDST 

SUBROUTINE HYST (SIGMA, GI) 

COMMON XI 

COMMON /CRK/ 00, XII, CG 
COMMON /CHS/ LAMBDA, PSO, PS 1,DX,RDX,CMA, V,W,PHI ,CMC,WR,RKEP 
COMMON /COR/ Ot GO 
REAL LAMBDA 

COMPLEX AA,AB,BA,BB,CA,CB,DX,GD,CG,TDX 
CG = PSO 

IF (SIGMA. NE.O.) CG = S IGMA*PS1* ( 0 . , 1 . ) 

C PSO = 0 WHEN BEARING IS WHIRLING 

C PS0=0 MEANS ORIFICE MASS FLOW DOES NOT VARY WITH THETA 
GD = ( 0 • , 1 ■ ) # ( LAMBDA - SIGMA) 

DX = • 5*X I 
TDX = 0. 

DO 110 N = 1, B 
BA = 0. 

BB = (-1 .,-1. ) 

CALL RKLP (BA,DX,AA,CA) 

CALL RKLP (BB,DX,AB,CB) 

CA = 1.57079633/(XI+XI1)*(AB*CA - AA*CB)/(AB - AA) 

IF (CABS( TDX/CA— 1. ) .LE. RKEP) GO TO 120 
DX = . 5 * DX 
110 TDX = CA 

120 RDX = REAL ( DX ) 

GR = —REAL ( C A ) 

GI = A I MAG ( C A ) 

W = CABS ( CA ) 

PHI = ATAN2(GI ,GR ) *57 .2957796 
WV = W/(V-1.) 

WR = 1000. 

CMC = 0. 

IF ( SIGMA. NE.O. ) CMC = 1 44 . * GR / S I GM A**2 
CMA = CMC*LAMBDA**2/36. 

CMS = CMA/ (V-l. ) 

IF (LAMBDA. NE.O.) WR = . 5*S I GM A/ L AMBDA 

WRITE (6,145) LAMBDA, GR, GI ,W,WV, PHI , CMA, RDX, CMC, WR, SIGMA 
145 FORMAT ( 3G 1 1 . 3 , 2 G 1 2 . 4 , G 1 1 . 3 , G 1 2 . 4 , 4G1 3 . 4 ) 

RETURN 

END 

$ I BFTC RKLPS 

SUBROUTINE RKLP ( G I N , DX , AD , C D ) 

COMMON XI 

COMMON /CRK/ QD, XII, CG 

COMPLEX G ( 4 ) , A ( 4 ) , B<4), XF, DX, AD, GIN, CD, CG 
EXTERNAL DLP,DLP1 
G ( 1 ) =0. 

G ( 2 ) - 0. 

G ( 3 ) = GIN 
G ( 4 ) = 0. 

XF = XI 


CALL RKGC ( DLP , G , XF, DX , A, B , A ) 

G ( 3 ) = G ( 3 ) + OD + CG*G<2) 

IF ( XI l.EQ.O. ) GO TO 10 
XF = XI + XII 

CALL RKGC ( DL P 1 , G , XF , DX , A , B , 4 ) 
18 AD = G ( 3 ) 

CD = G ( A ) 

RETURN 

END 

SIBFTC DRVLP 

SUBROUTINE DLP(G,GP) 

COMMON /C DR/ Q, GD 
COMMON XI 

COMPLEX GD » G ( 4 ) ,GP(4) 

T = REAL (G( 1 ) ) 

GO TO 10 

ENTRY DLP1 (G,GP) 

T = XI 

10 GP ( 1 ) = 1 . 

GPU) = G ( 3 ) 

S = SORT ( 1 ,+0*T ) 

GP ( 4 ) = G(2)/S 

GP( 3) = G( 2 ) - GD* ( GP ( 4 ) + S) 

RETURN 

END 


SIBFTC RKGHC 

SUBROUTINE RKGC (DERIV, Y, XF INAL * DELTA, 0, YP, N) 

C RUNGE-KUTTA-GILL INTEGRATION OF N-l COMPLEX FUNCTIONS OF COMPLEX 
C ARGUMENT Y(l) FROM Y(l) INITIAL TO XFINAL 
DIMENSION A ( A ) , B ( A ) , C(A) 

COMPLEX Y(l), 0(1), YP(1), XFINAL, DELTA, T 

DATA A/,5, ,292893219, 1.70710678, .166666667/, B/2., 1., 1., 2./, 

1 C/.5, .292893219, 1.70710678, .5/ 

KK = CABS( (XFINAL -Y(l) )/DELTA) + .01 
IF (KK.LE.O) GO TO 20 
DO 5 I = l v N 
5 0(I)=0. 

DO 10 K = If KK 
DO 10 J = 1, A 
CALL DERIV(Y,YP) 

DO 10 1 = 1, N 

T=A( J )*( YP( I )-B( J)*0( I ) ) 

Y( I )=Y( I ) +DELTA*T 
10 Q(I)=0( I )+3.*T-C(J)*YPI I ) 

20 RETURN 

END 
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P— PERTURBATION 

ANALYSIS FOR 

HYBR10 BEARING IN STEADY 

CIRCULAR 

WHIRL 

L I/O 

LAMtiDA T 

1+DELTA**2 

K=CP/CV 

ALFA 

V = PS/ P A 

l 0.750 

4.000 

1.200 

1.400 

0.700 

5.000 


SOLVE FOR YA = PC/V 

1 YA 0.528282 DYA 0.459088E-01 

2 YA 0.574191 DYA - 0. 2 85489E- 02 

3 YA 0.571336 OYA -O. 1 12002E-04 

4 YA 0.571324 DYA 0.I86657E-08 


PC/PS 

0.571 

PSO 0.640E-01 

PSl 0. 

210 RKOX 

REOJCEO UNTIL ERROR 

.LE. 

0.500E- 

02 




LAMBDA 

FR/EPALO 

FT/EPALO 

M/EPALD 

W/E ( P S-PA ) LD 

PHI MCW2/PALD 


R< DX 

MPC5/2LI5MU2 

WHIRL RATIO 

VI 8 

NR SI3MA 

0 

2.352 

-0.313E-06 

2.3521 

0.5880 

-0.762E-D5 

0 


0. 1875 

0 

1000.0 


0 

0 

2. 389 

— 0. 1 87E-0 1 

2. 3890 

0.5973 

-0.450 

0 


0. 1875 

38.223 

1000. 0 


3.0000 

0 

2.415 

0.108E-02 

2.4149 

0.6037 

0.256E-01 

0 


0. 1875 

347.74 

l 000.0 


1.0000 

0 

2.414 

0.856E-03 

2.4138 

0.6034 

0.203E-01 

0 


0. 1875 

282.73 

1000.0 


1.1088 

0 

2.409 

-0.909E-03 

2.4088 

0.6022 

-0.216E-Q1 

0 


0.1875 

148.44 

1000.0 


1.52S6 

0 

2.411 

0. 187E-03 

2.4114 

0.6029 

0 . 444E-02 

0 


0. 1875 

201.60 

1000. 0 


1.3124 

LAMBDA 

FR/EPALD 

FT/EPALD 

W/EPALD 

W/ECPS-PA1 LD 

PHI MCW2/PALD 


RKD< 

MPC5/2LR5MU2 

WHIRL RATIO 

VI B 

NR SIGHA 

0. 500 

2.362 

0.116 

2. 3647 

0.5912 

2.B15 

0 


0.1875 

0 

0 


0 

0.500 

2.336 

0.623E-01 

2. 3365 

0.5841 

1.529 

0.2595 


0. 1875 

37.370 

3.0000 


3.0000 

0.500 

2. 313 

—0. 1 1BE-Q1 

2. 3125 

0.5781 

-0.292 

0.9250E- 

01 

0. 1875 

13.320 

5.0000 


5.0000 

0.500 

2.314 

-0.160E-03 

2. 3140 

0.5785 

-0.397E-D2 

0.1056 


0. 1875 

15.199 

4.6822 


4.6822 

0.500 

2.314 

0. 159E-05 

2.3140 

0.5785 

0. 3936-04 

0.1058 


0. 1875 

15.228 

4.6778 


4.5779 

LAMBDA 

FR/EPALO 

FT/EPALD 

W/EPALO 

W/E1PS-PA1L0 

PHI HCW2/PAL0 


RXDX 

HPC5/2LR5MU2 

WHIRL RATIO 

VI B 

NR SIGHA 

1.000 

2.390 

0.228 

2.4010 

0.6003 

5.439 

0 


0. 1875 

0 

0 


0 

1.000 

2.289 

0.156 

2.2943 

0.5736 

3.902 

1.0173 


0. 1875 

36.624 

1.5000 


3.0000 

1.000 

2. 252 

0.466E-01 

2.2530 

0.5632 

1. L85 

0.3604 


0.1875 

12.974 

2.5000 


5.0000 

1.000 

2.258 

-0.443E-01 

2.2587 

0.5647 

-1.L25 

0.1844 
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Figure 8. - Example of working sheet 
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Figure 9. - Calculation summary sheet 
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